home *** CD-ROM | disk | FTP | other *** search
/ MacFormat 1995 January / macformat-020.iso / Shareware City / Education / RLaB / examples / sfft.r < prev    next >
Encoding:
Text File  |  1994-02-21  |  357 b   |  26 lines  |  [TEXT/ttxt]
  1. //
  2. // Calculate Fourier transform of a finite signal (SLOWLY)
  3. //
  4.  
  5. sfft = function ( x )
  6. {
  7.   local(i, k, n, X);
  8.  
  9.   if(class (x) != "num") 
  10.   {
  11.     error ("input to sfft must be a matrix");
  12.   }
  13.  
  14.   X = zeros ( size (x) );
  15.  
  16.   for( k in 1:x.n )
  17.   {
  18.     for( n in 1:x.n )
  19.     {
  20.       X[k] = X[k] + x[n]*exp( (-1i*2*pi*(k-1)*(n-1))/x.n );
  21.     }
  22.   }
  23.  
  24.   return X;
  25. }
  26.